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Abstract 

I use a previously introduced mapping between the continuum per- 
colation model and the Potts fluid to derive a mean field theory of 
continuum percolation systems. This is done by introducing a new 
variational principle, the basis of which has to be taken, for now, as 
heuristic. The critical exponents obtained are (3 = 1, 7 = 1 and 
v = 0.5, which are identical with the mean field exponents of lattice 
percolation. The critical density in this approximation is p c = 1/VJ 
where VJ = J dxp(x){exp[— v(x)/kT] — 1}. p(x) is the binding prob- 
ability of two particles separated by x and v(x) is their interaction 
potential. 
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1 Introduction 



The first paper in this series [l|], hereafter referred to as I, contained a general 
formalism of continuum percolation. Such a system consists of classical par- 
ticles interacting through a pair potential v(fi,fj), such that they can also 
bind (or connect) to each other with a probability p(fi,fj). The function 
qifi^Tj) = 1 — p(fi,fj) is the probability of disconnection . 

The formalism is based on a quantitative mapping between the percola- 
tion model and an extension of the Potts model I have named the Potts fluid. 
For easy reference, I recall here the essential definitions and results. 

The s-state Potts fluid is a system of N classical spins {Aj}^ interacting 
with each other through a spin-dependent pair potential V(fi, A^; r}, Aj), such 
that 

V(f„K;f 1 ,X i )^V ( i, j )^[^ ) !j £ = £ (1.1) 

The spins are coupled to an external field h(r) through an interaction Hamil- 
tonian 

N 

H int = -Y,^i)Wi) (1-2) 

i=i 

where 

Up to some unimportant constants, the Potts fluid partition function 
(more precisely, the configuration integral) is 



Z = Jtf\ H / drf- df N exp 



{A m } 

The magnetization of the Potts fluid is defined as 



N 

0£v(vj) + 0EMOV>(Ai 

i>j i=l 



;i-4) 



1 9\nZ 

where h is the now constant external field. The susceptibility is 

dM (i a\ 

X = ^h (L6) 
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The n-density functions of the Potts fluid, are defined as 
p (n) (fi,Ai;f 2 ,A 2 ;...,f n ,A n ) = 

N 

-/?5>( < ».7')-/?EMiMAi) 



- / dr n+ i ■ ■ ■ df N exp 

— n ! J 



Z(N -n)\ 



i>j i=l 



(1.7) 

Of particular interest is the spin pair-correlation function, defined as 

gi 2) {x,a;y,-f) = 1 p (2) (f , a; y, 7) (1.8) 
p{x)p{y) 

which tends to 1 when \x — y\ — > 00. Here p(x) is the numerical local density 
at the point x. It is often useful to define a connected spin pair-correlation 

as 

hf\x,a;y,7) = gf ] '(£,«;?/, 7) - 1 (1.9) 

This function tends to zero when \x — y\ — > 00. 

Any continuum percolation model defined by v(i,j) and p(i,j) can be 
mapped onto an appropriate Potts fluid model with a pair-spin interaction 
defined by 

U(i,j) = v(i,j) 
exp[-PW(i,j)] = q(i,j)exp[-(3v(i,j)} (1.10) 

The relation between the Potts magnetization and the percolation prob- 
ability, -P(p), is 

lim lim lim M = Pip) (1.11) 

h— >0 N^oo s— >1 

For densities lower than the critical density, the susceptibility is directly 
related to the mean cluster size, S. 

lim lim lim X = PS (p < p c ) (1.12) 

/»— >0 >oo s— >1 

An important quantity in the percolation model is the pair-connectedness 
function g'(x,y), the meaning of which is 

p(x) p(y) g\x \ y) dxdy = Probability of finding two particles in regions 

dx and dy around the positions x and y, such 
that they both belong to the same cluster. 

(1.13) 
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This function is related to the mean cluster size by 



S = 1 + p J dfg\r) 



;i.i4) 



where we assume, as we usually shall, that the system is translationally 
invariant, so that g^(x, y) = g\x — y) and p(x) = p{y) = p. 

The pair-connectedness is related to the Potts pair-correlation functions 

by 



lim 
lim 



gf\x, o\ V, cr) - 9i 2) (x, a; y, rj) 
hf ] (x, a; y, a) - hf ] (x, a; y, rj) 



;i.!5) 



where a, i] ^ 1 and o ^ r\. 

As usual, the first approach to any phase transition is to find the mean 
field theory. This is the aim of the present paper. The mean field theory for 
the magnetization (i. e., the percolation probability) and for the susceptibil- 
ity (z. e., the mean cluster size) is developed in section 2 from a physical point 
of view. Then, in section 3, I show that the same results can be obtained 
from a variational principle, which however, remains at present a heuristic de- 
vice. Section 4 uses the variational principle to obtain the pair-connectedness 
within the mean field approximation. Section 5 sums up the results. 



2 Mean field theory for M and % 

As usual, the basic assumption of mean field theory is that every spin A feels 
an average, homogeneous interaction H (A) / (3 due to all the other spins. Since 
the field h distinguishes spins in state 1 from all the others, we assume that 
H(\) takes two values, if (A = 1), and if (A ^ 1), this last being identical for 
all spin states other than 1. This is because the field preserves the symmetry 
between the non-1 spin states. 

From Eq. (|1.4j) , the configuration integral (or the partition function, up 
to nonimportant constants) is, in this approximation, 

if ) N 
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definition of ^(A), Eq. ( |1.3|) , we have 



Separating now the case A = 1 from the s — 1 cases A 7^ 1, and recalling the 

(2.2) 



dfn 



N 



with 



Q = exp[-H(\ = l) + /3(s-l)h(r)] 
+ (s - 1) exp [-PYA ^ 1) - /?/i(r)] 



Let us denote 



B v = exp [-H(X = 1)] 
B r = exp \-H(\ 1)] 



For a constant field, Eq. ( j2.3|) becomes 

yN 



N 



Bie P(s-i)h + ( s _ i) Br6 -Ph 

The magnetization is obtained from Eq. ( |1.5|) , 

1 d\nZ Btfft'- 1 * - B r e~P h 

M 



N(s - l)/3 dh B^- 1 ^ + {s- l)B r e~P h 
In the percolation limit, s — > 1, we have 



M = 1 — lim 



5, 



or, for h = 0, when M becomes the percolation probability P, 



P = 1 - lim 



P, 



1 VPi 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



This equation is a self-consistency condition because (B r /Bi) depends on 
M. The reason is as follows. A spin 1 interacts differently with other spins 1 
than it does with spins in other states. As a result, Pi must depend on the 
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average number of spins in state 1 and of spins in all other states. The same 
reasoning holds for B r . Now, let us denote 

ni = number of spins in state 1 

n = number of spins in any state a/1 (2-9) 
I have shown in I [Eq. (4.9)] that 

M=^{m-n) (2.10) 

Therefore B\ and B r depend on M, (or P) through n\ and n. 

Finding the form of (B r /Bi) requires several approximate arguments, 
which will be clearer if we start with a simple case, i. e., v(i,j) = 0. In 
the percolation picture, this means no interactions . In the Potts fluid pic- 
ture, this means that the only interactions take place between non-parallel 
spins, and these interactions are given by exp[—/3W(i,j)] = q(i,j) [compare 
Eq. ( |1.10| )1. In this case, a spin 1 interacts only with the (s — l)n spins 
which are in non-1 states. Similarly, a non-1 spin interacts only with the 
ni + (s — 2)n spins which are non-parallel to it. 

The first step is to notice that here we can take already the limit s — > 1, 
since the term (s — 1) in the denominator in the definition of M, Eq. ( |1 . 5| ) , 
has vanished in Eq. ( |2.7| ). In this limit, a spin 1 interacts with zero [the limit 
of (s — l)n] other spins, while any non-1 spin interacts with rii — n = NM 
other spins. 

The main argument is now that it is easier to find B r and B\ in the limit 
M — > 1, i. e., when practically all the spins are in the state 1. 

Consider first a state where all the spins are in the state 1. Such a state 
can be thought of, for example, as the limit of an infinitely strong field. Then 
there are no inter-spin interaction, and the partition function is 

1 r V N 
Z = — J dl---dNexp[pNh(s-l)} = —exp\pNh(s-l)} (2.11) 

However, by definition, we also have 

yN 

Z=—B?exp[(3Nh(s-l)] (2.12) 

Comparing the equations, we have that B\ = 1. In principle, this holds in the 
limit M — > 1 only; however, as mentioned above, in the s — > 1 limit, a spin 
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1 interacts with no other non-1 spins. On the other hand, in our case, the 
interaction with other spins 1 is v(i,j) = 0. Therefore, all the interactions 
vanish and the above result for B\ can be extended to all values of M. Hence, 
B\ = 1 in general. 

Next consider a state where all spins but one are in the state 1 (in the 
thermodynamic limit, this is still M — 1 ). Again, the partition function can 
be calculated easily. For defmiteness, assume that the non-1 spin is always 
numbered N. Then, from Eq. (|1.4j) , we have 



1 



N-l 



(N-ly.I dl "' dN ^J-^E^Ml + W- 1 ) 5 -^ 
i 



n N-l 



(N-l)! 
And we must also have 

yN 



dnq(f N - n) 



exp {(3h[(N -l)s- N]} 



(N-l)\ 



BrB?- 1 exp {(3h[(N - l)(s - I) - 1]} 



Hence, 



B, 



— J drq(r) 



N-l 



Using p(f) = 1 — q{f) and denoting 



we have 



B r 



dfp(r) , 



1 - 



V 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



Now the power N — 1 represents the number of spins which interact with the 
non-1 spin. In other words, it is n\ — n for this particular configuration. In 
the thermodynamic limit it can therefore be replaced by MN. In terms of 
the density p = N/V we can write 



B, 



N 



-, MN 



(2.18) 



Now although this equation was derived for the case M — 1, it is actually 
written for an arbitrary M. Consequently, we'll assume now that we can use 
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it for all values of M. This is part of the mean field approximation itself, 
i. e., that for the case of v(i,j) = 0, we suppose that in the thermodynamic 
limit, 

pV e ^ MN 



— = hm 

B 1 N->oo [ N 

Substituting this into Eq. (|2.7|) yields 



exp (-MpV e ) (2.19) 



M = 1 — exp {-MpV e - fish) (2.20) 

or, for h = 0, 

M = 1 — exp (-MpK) (2.21) 

This equation is very similar to the one obtained in the mean field theory of 
the Ising model. For low densities, its only solution is M = 0. A non-zero 
solution appears first when the slopes of the two sides of the equation are 
equal at M = 0. This condition yields for the critical density, p c , the result 

Pc = y £ (2.22) 

Typically 0, the function p(f) is chosen to be 

1 r belongs to some excluded volume, V exc , 
p(r) = ^ around the center of the particle (2.23) 

else 

Then, V e = V exc , and Eq. ( |2.22| ) becomes 

Pc = tt~ (2.24) 

'exc 

Computer simulations show this result to be correct in the limit of infinite 
spatial dimension. This is precisely where mean field theory is expected to 
be correct. 

Expanding now Eq. (|2.2U ) in the vicinity of the critical point, we have 
M = M-^-^M 2 +... (2.25) 
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Hence, the percolation probability, which is equal to M in this limit, is 



P{p) = M 



2p c 



(Pc ~P)~ {Pc - pf {p -> Pc) 



(2.26) 



where in this equation, (3 is the critical exponent of P(p). Therefore, we find 

(3 = 1 (Mean Field) (2.27) 
We can now calculate the susceptibility. From Eq. (|2.20|) , 

dM 



X 



dh 















h=0 


h=0 



exp {-MpV e 



or 



X 



0e 



-MpVe 



(2.28) 



(2.29) 



l-pV e 

For p < p c , we have that \ — > f3S. In this case, M = 0, so that the mean 
cluster size is 



X 



S 



1 



P l-pV e 
Hence, the critical exponent 7 is 



{pc - py 



7 = 1 (Mean Field) 



(2.30) 



(2.31) 



We now wish to extend this argument to cases where v(i,j) 7^ 0. The 
essentials remain unchanged. Again we look at the limit M — > 1. This time, 
when all the spins are in the state 1, the equivalents of Eqs. (|2. lip and ( [2.12j ) 
yield 

1 r T iv 

£?f = — Jdl---dNex P 



i=l i>j 



(2.32) 



Similarly, if only one spin is in a state other than 1, we have 
B^Br = — J dl ■ ■ ■ dN exp 



AT-l 
i=l i>j 



N-l 



x U q(i,N)exp[-(3v(i,N)) (2.33) 
i=i 
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Dividing the two equations, we find that in this limit 



i. e. 



M -> 1) 



B r 



N-l 


N-l 


Jdl---dN n q{i,N)exp 


-P E v(i,j) 


-Pv(i,N) 


i=l 


i>j=i 




N 




f dl ■ ■ ■ dN exp 


-13 E «M 






i>i=i 





(2.34) 



In general, one cannot compute exactly this ratio. However, we don't need 
B r /Bi exactly, but only in the region M ~ 0, where the critical behavior 
occurs. Now, from the previous arguments, B r /Bi can depend on M only 
through the excess spin-1 density Mp = p 1 — p a (where p 1 is the density of 
spins 1 and p a is the density of spins in any state a / 1). Hence, we are 
interested in the region of small Mp. At first sight, this seems unhelpful, 
because Eq. (|2.34j ) corresponds to the limit M — *■ 1. However, Mp may be 
made small by decreasing the density rather than M. In other words, the 
region of interest Mp — ► may be obtained by looking at the limit p — » 
even if M — > 1. 

We therefore take the small density limit of Eq. ( |2.34j ) and extract the 
leading behavior. Let us introduce the functions 



f(r) = exp[— /3v(r)] — 1 
f*(r) = q(f) exp[— (3v(r)] — 1 
f\r) = /(f) - f*(r) = p(r) exp[-/3v(r)} 



(2.35) 



The function f{f) is just the Mayer f-function ||. Its usefulness stems from 
its being typically short ranged (i. e., f(r) — > at |f| — > oo). Therefore, inte- 
grals over / (and over /* and f> as well) remain finite in the thermodynamic 
limit . 

Eq. ( |2.34|) now becomes 



B r 
B~i 



N-l 

Jdl---dN n n [l + f(iJ)][l + f*(i,N)] 

i=l [*>j 

N-l 

jdi---dN n n [i + f(hj)\ 

i=l i>j 



(2.36) 
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Expanding this expression, we find 



B~i 



Jdl---dN< 


N-l 

i 1+ £ 




*■■■} 


fdl- 


• ■ dN < 


AT-l 

1+ E E/(*,i) + "- 

2=1 t>j 


1 



1 + [(JV- l)(iV-2)/2V]/df/(f) + (AT — 1/V) / dff*(r) + 
1 + [(AT — 1)(AT — 2)/2^] / drfif) + (AT — 1/V) J dff(f) + ■ 



(2.37) 



Keeping now only the leading terms in the density, we obtain 



Br „ 1 + {N- l/V)Jdrf*{f) + 
~B~ X ~ 1 + (N- l/V)Jdff(f) + • 



1 - 



r 



drf{r) + --- (2.38) 



where we have used the definition f'(r) = /(f) — /*(f)- Again, as in 
Eq. ( |2 . 1 8|) , we recognize that the term AT — 1 represents just n\ — n in this 
configuration. Hence once again, we can replace it with MN and write 

MN r + + 
' df p(f) + ■ ■ ■ = 1 - MpVj 



V 



where 



v; f = / df / f (f) 



generalizes Eq. ( |2.40|) , to which it properly reduces if v(r) = 0. 
To the same order in p, we may also express Eq. fl2.39|) as 

B r 
B[ 



exp 



-MpV} 



(2.39) 
(2.40) 

(2.41) 



which completes the analogy with Eq. fl2.19Q for the case v = 0. Because of 
this, all the results for the case v = apply to the more general well. 
Hence, for arbitrary potentials and binding criteria, we now find 

P = 1 - exp (-pPV}) 

(P < Pc) 



s 

Pc 
IMF 



1 - pVj 

1 _ 1 

v?~ WW) 
i 

i 



(2.42) 
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3 A heuristic variational formulation 



We would like now to calculate the connectedness function in the mean field 
approximation. The previous arguments are inadequate for this task. One 
knows, however, that in the Ising model, the correlation function may be 
calculated within mean field theory through a variational formulation based 
on the Bogolyubov-Feynman inequality H, i. e., 



> exp (-(3(Hi)o) 



(3.1) 



where the subscript refers to some reference system and Hi = H — Ho, 
where H is the true Hamiltonian of the system and H® is the Hamiltonian of 
the reference system. 

Unfortunately, this inequality is inadequate to deal with our system, be- 
cause the potential v(i,j) will typically have a strongly repulsive part ("hard 
core") at short distances, or else the function q(i,j) will typically vanish in 
some range (where the binding is certain), both of which make the r.h.s of the 
inequality undetermined. One would therefore like a variational formulation 
of the mean field which allows for strong effective interactions. I will show 
that the mean field developed in the previous section can be derived from 
such a principle, though I have been unable to formally prove the principle 
itself. Its status must therefore be taken at present to be essentially heuristic. 

The motivation is as follows. Let us assume, as usual, some reference sys- 
tem in which the spins interact only with an external field B. The reference 
Hamiltonian, Ho, is 

JV 

H = -Bj2?P(k) (3-2) 
i=i 

Then the partition function can be written as 



— J2 I cFi • • • df jv | exp 
iV '{A m } J I 



2V 



8=1 



Z ( exp 



x exp 



i>j i=l 
N 



N 



i=l 



(3.3) 
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where 



Z o = Tn E [ dn - ■ ■ dr N exp 

iV ' {A m } J 



iV 



i=i 



A' 



V 
AH 



+ (s-l)e 



A? 



(3.4) 



and ( )o means an average performed in the reference system with the 
Hamiltonian Tio- The convexity of the exponential function then gives us 
the Bogolyubov-Feynman inequality, 



> exp 



N 



-{3j2V{i,j) + P{h- B^iPiXt 



i>3 



i=l 



(3.5) 



The usual mean field theory follows from extremization with respect to the 
parameter B. The resulting approximation to Z is expected to be good as 
long as (3V is small. Now, if this is the case, we could as well write 



-pV(i,j)nex V [-/3V(i,j)]-l 



so that to this order 



— ~ exp 



~/3V(i,j) 



N 

l) Q + i3(h-B)J2mi))o 

1=1 



(3.6) 



(3.7) 



However, at this stage it is not clear what has become of the inequality 
in Eq. (|3.5| ), because the replacement Eq. (|3.6|) is increasing the exponent, 
thereby countering the Bogolyubov-Feynman inequality. 

Nevertheless, for small /3V, we expect that (l/N)lnZ can be approxi- 
mated by the following function, 



1 



N 



F t (B) = - In Z + £ (exp [~pv{h j)\ ~ 1>„ + E P( h ~ B ) &( X i))a 
iV [ i>j i=i 

(3.8) 

Typically, of course, V(i,j) is not uniformly small. The surprising result, 
however, is that if we use for (l/N)lnZ the extremum value of F t (B) with 
respect to B, we shall recover the mean field theory of the previous section. 
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The proof is straightforward. First, 
WA .)) EA^(Ai) exp [/3^(A,)] 



1) 



3 /?B(s-l) 



e /3B(s-l) + ( s _ l) e -/9B 



E Ai exp [/3S^(Ai)] 
Next, we rewrite 

ex P [—0V(i, Xf, j, Xj)] - 1 = j)exp[-/?f(i, j)] - 1 



(3.9) 



(3.10) 



Taking the average with respect to the reference system of the two terms on 
the r.h.s, we find first 

(q(i,j)exp[-f3v(i,j)} - 1) = ^ J didj {q(i, j) exp[-(3v{i, j)] - 1} (3.11) 

which is independent of B. 
The second term is 

p(i,j)exp[-f3v(i,j)]5 Xi ,\ j 



— [ dl---dN P(hj) ex V 
oI ^ lJ {x m } 



N 



-{3v(i,j)+PBj2^(Xk 



k=l 



^ I'lidj !>(>■.))<- Ir "" 



E 6 Xi>x .e W {pB[i(){X i )+il){X j )}} 



#B{a-l) + ( s _ iy 



V 2 



didj pfafie-WA 



- l)e- 2 ? B 



^B(s-l) + ( s _ !) ( 



Note that from the definition of V^, Eq. (|2.16 ), we have 



V 2 

Substituting Eqs. Q3.4Q , 
Eq. (|3.8|), we have 



didj p(i, j)exp[-/3w(i,j)] 



(3.12) 



(3.13) 



and ( |3.12|) into the definition of F t (B), 



Ft(B) 



C + \n 



e?B(.-i) + r s 



l)e 



-I3B 



+ 



iV- 1 /y e t\ e 2 ^" 1 ) + (s- l)e- 



2/3B 



V 



+ (3(h-B)(s-l) 



V J [ef>B{a-l) + ( s _ l) e -/9B] 
e /8-B(s-l) _ e -/9B 



e^B(»-i) + ( s _ i) e - 



/3B 



(3.14) 
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where 



C = (q(i,j)ex V [-Pv(i,j)\ - 1) + \nV - - ln(iV!) (3.15) 



is independent of 5. 

We now have immediately that 



dFt _ (3 2 s 2 (s- l)(h-B)e 
~dB ~ 



-/3Bs 



[1 + (s-l)e- 



-PBs 



+ 



(n-i)v: 



V 



Ps(s - l)e~P Bs (l - e~P Ba ) 

[1 + (s - l)e-? Bs f 

(3.16) 



The extremum condition dF t /dB = yields the equation 



B — h 



'pvy 



i 



-PBs 



[s - l)e-P Bs 



(3.17) 



In this approximation, the magnetization is 
1 dF t 1 



M 



(3(s-l)dh P(s-1) 



dF t OB + (dF t 



OB dh 



dh 



'dF t 



(3(s-l) \ dh 



where we have used the extremum condition dF t /dB = 0. 
From Eqs. (|3.14|) and ( |3.17 ), we have 



M 



eflB<a-l) _ e ~PB f3s 



e PB(s-i) + ( s _ \) e -pB pVe 
Hence we can rewrite Eq. (|3.17| ) as 

(3sB = fish + MpVj 
Substituting this result back into Eq. ( |3.19| ) gives us 

I _ e -M P vJ e -p s h 



j(B-h) 



M 



1 + (s - l)e- 



-MpVj e -Psh 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



which is exactly the mean field equation, Eq. ( |2.6p with B r jB\ given by 
Eq. ( |2.41D . Hence, indeed, all the mean field results may be obtained from 
the variational principle of F t . 

Finally, note that from Eq. (|3.16 ), the extremum of F t is actually a max- 
imum. This might suggest the existence of an inequality InZ > NF t , but I 
have been unable to prove it one way or another. 
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4 Correlation functions 



The variational principle now allows us to find the correlation function in 
mean field theory. To this end, we need to relate the Potts n-density functions 
defined in Eq. (|1.7| ) to the partition function. By analogy with a similar 
formalism in the theory of liquids ||], let us define a generalized functional 
differentiation operator, in the following way. Let T\t{f, A)] be a functional 
of the function t(f, A), which depends on a position variable as well as on an 
associated discrete spin variable. Then the generalized functional derivative 
SJ-'/dt will be defined through the relation 



j^^WA) (4.1) 



where 8T is the change in T associated with a variation 5t in t(r, A). The 
only difference with the usual functional derivative is in the added summation 
over the spin variable. This does not change any of the basic properties of 
the operator. 

I shall show later that the Potts correlation functions are functional 
derivatives of the partition function, for fields h which are inhomogeneous. 
In order to extend the mean field theory to cover this case, we now need to 
adopt the following heuristic claim: 

Heuristic Assumption: The F t variational principle yields the proper 
mean field theory of continuum percolation even for the case of inhomoge- 
neous fields h(f) and B(r). The extremum condition must then be general- 
ized to be 

*3. = U 2) 

which determines the extremizing field B(f). 

The application of this extremum condition yields a complicated integral 
equation. However, for our purpose, we only need F t [B] in the region h{f) — > 
and in the vicinity of the critical region M(r) — > 0. Looking at the results 
for the homogeneous case, Eq. (|3.20|) , we can expect that the extremizing 
function B(f) will be also very small in this case (the solution shows this 
assumption to be self-consistent). Therefore, one only needs to look at the 
extremum condition for very small B(r). The details of this calculation are 
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in the Appendix. The final result is the condition (to first order in B) 

B(x) = h(x) + ^Jdyp{x- y)e-^ 3 ~^B{y) (4.3) 

where we have assumed that p(x, y) = p(x — y) . Note that if B and h are 
uniform, we recover the first order limit of Eq. fl3.17| ), as we should. This 
equation is easily solved in Fourier space where 



B(k) 



h(k) 



(p/s)P(k) 



(4.4) 



with 



f\f) = p(r)e~^ (4.5) 

and functions with hats are Fourier transforms of the corresponding functions 
without hats. 

Let us now calculate the Potts correlation functions. We take 



t(f,X) = exp[/3h(r)ip(X)] 
Then we have for the partition function from Eq. (IT 



1 r N 

z = w £ / ^ • • • df N n y(ij) 



{A m } 

It is now simple to see that 
~5Z N 



i=i 



N 



df 2 --- df N Y[ tin, Xi) 



i=2 



Comparing this to Eq. ( |1.7|) we see that 



P {1) (r,A) 



t(f, A) SZ 



5\nZ 



Z 5t(f,X) 6hit(f,\) 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



Similarly, we have that 
8®Z 



6t(x,a)6t(y,7) 
1 



N 



df 3 --- df N Yl t(n, Xi 



i=3 



r2=y; x 2 =i 



(4.10) 
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And a comparison with Eq. ( |1 . 7| ) shows that 

P (2) (x, a; y, 7) - p (1) (x, a)p {1) (y, 7) = t(x , a)t(y, 7) — , , , 

dt(x,a)dt(y,j) 



<5 2 lnZ 



(4.11) 



Within the mean field theory we replace InZ with NF t [B]. Making use 
of the extremum condition, we now have from Eq. Q4.9|) and the definition of 
t(f, A), Eq. flU), that 



P W (x,a) 



SNF f 



5[(3h{x)^(a)} 



N 



exp [pB(x)ip(a)] (4.12) 



where 



r = [dzexp [pB(z)ip(\)\ 



(4.13) 



The easiest way to calculate the pair-correlation function is now to use the 
identity 



Sp^'(x, a) ._ „ 5 

= t{y n ) 



5[(3h(y)^)} 



= t{y,i) 



St(y,j) 
5NF t 



t(x, a) 



5NF t 
5t(x, a) 



S(x-y) S an + t(x , a)t(y, 7)- 



5 2 NF t 



'5t(x,a) v 1 ' §f^x, a)8t(y, 7) 

= (£, - y) 5 a ^ + p (1) (£, a)p (1) (y, 7)/i (2) (£, a; y, 7) 

(4.14) 

On the other hand, we have from Eq. ( |4.13| ) that 



6pV>(x, 



a) 



5[(3h(y)^)} 



1 - -e 



/3B(j?)i/>( 7 ) 



r mmi)} 1 ' ] 



We are interested in this function in the range p < p c and h = 0, which 
is the one relevant for percolation. In this case, both h and B vanish, so 
that we can evaluate all quantities at zero fields. Thus, using the fact that 
T = sV in this limit (see the Appendix), we have that 



p^\x,a] 



_ N 

b=o r 



= p 

B=0 S 



(4.16) 
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where p = N/V as usual. In the same limit, we also have from Eq. ( [4.15| ) 

(4.17) 



Sp^(x, a) 


_ P 




5[B{x)ip{a)} 






h=0,B=0 S 


sV. 




h=0,B=0 



In the thermodynamic limit, the second term inside the brackets on the r.h.s 
vanishes (its presence is typical of the use of the canonical rather than the 
grand canonical ensemble). Combining Eqs. (|4.17f ), Q4.16Q and ( 4.14Q , we 
obtain 



P 4h^x,a;y^)^ P -~ 5[B{ma)] 



s 6[h(y)4>(7)} 



h=0,B=0 



- v) 8, 



Q,7 



(4.18) 



Since the system is now homogeneous, hP' depends only on the difference 
x — y . Taking the Fourier transform of Eq. ( [4.1 8|) , we have that 



Pr 

--o, 



a, 7 



- I dke 



ik(3-y) S[B(x)jj(a)} 
5[h(y)^)} 



(4.19) 



h=0,B=0 

Finally, we note that from the definition of the generalized functional deriva- 
tive, 

5[h(mi)\ &KD ' 7 1 j 

where the expression 5B/5h on the r.h.s is now a usual functional deriva- 
tive. The properties of the usual functional derivative imply that 



SB(x) 



5h(y) 



5B(k) 



h=0,B=0 



Sh(k) 



h=0,B=0 



(p/s)P(k) 



(4.21) 



where in the last equality we have used Eq. (|4.4|). Substituting this result in 
Eq. Q4T9D finally yields 



/i( 2 )(£;a, 7 ) 



-5, 



a, 7 



(4.22) 



1 - (p/s)p(k) 

Hence, substituting this into Eq. ( |1.15|) finally yields the pair connectedness 
within the mean field approximation, 



i- P p(k) 



(4.23) 
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This expression is consistent with Eq. ( |2.42| ) for the mean cluster size through 
the relation between S and given by Eq. (|1.14|) . 

Since f'(r) = f '(— r), then for small k, we have that 

/t(jfe) = Vl - Vlk 2 + 0(k 4 ) (4.24) 

where 

V} = J dff\r) 

2 

V? = jdr r -f\r) (4.25) 
Therefore, for small wavevectors, the pair-connectedness is 



~ ' i - pv e j ) i + \y}/{\ - P v2)\w (4 ' 26) 



We expect by analogy with other critical phenomena that 

*W ~ ttW (4 ' 27) 

where ^ is the connectedness length, which plays in percolation theory a role 
analogous the the one the correlation length plays in other critical phenomena 
IJ. Hence, we have obtained an expression for the connectedness length 
within mean field theory, 



1 - pVl Pc- P 

Thus, we find that the mean field value of the exponent v is 



e = ~ ~^T7f = ~ (Pc - pY V (4-28) 



v = - (Mean Field) (4.29) 

5 Discussion 

The results of the mean field theory are summarized below. 

P = l-exp(-pPV^ 
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1-pVJ 



t 



1 



1 




Jdrp(r) 



PMF — 1 



1 



IMF = 1 
1 

v mf — - 



(5.1) 



The critical exponents obtained are identical to those obtained in perco- 
lation on a Bethe lattice, usually considered as the mean field equivalent for 
the lattice percolation model. This therefore addresses one of the questions 
presented in the introduction of I, namely the extent and nature of the contin- 
uum percolation universality class. As far as mean field theory is concerned, 
we have obtained theoretical proof that this universality class encompasses 
all interactions v(i,j) and all binding criteria p{i,j) and that it is identical 
to the universality class of lattice percolation. This last, incidentally, can be 
considered as a particular type of continuum percolation with interactions 
fixing the particles on the sites of a lattice. Therefore, the question of the 
independence of the critical exponents on the interaction is connected essen- 
tially to the relation between continuum and lattice percolation. We know, 
however, that mean field theory always grossly overestimates universality, as 
demonstrated by its being insensitive even to spatial dimensionality. Thus 
one cannot give too much weight to these conclusions. Nevertheless, the 
mean field result is an encouraging first step. Certainly had we found any 
difference between continuum and lattice percolation in the mean field ap- 
proximation, or any dependence on the interaction or the binding criterion, it 
would have been highly improbable that a more elaborate calculation would 
have restored universality. Thus at least, our results are encouraging, if far 
from definitive. 

As to the critical density, we do not expect it to be quantitatively ad- 
equate, as mean field results never are. However, as the dimensionality of 
the system increases, these results become better and better. Thus, Alon 
et al. || have shown that indeed, for a system of permeable hyper-cubes, 
pc^exc —>■ 1 for high dimensions. 
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As another example, for a /^-dimensional system of hard-cores and soft- 
shells |7|, §, we have 

r(l - J) < n otherwise (5 ' 2) 

^')H ~ otherwisT'' < d (5 ' 3) 

where a and d are two length parameters. Then we obtain in the mean field 
approximation 

1 _ 1 , . 

Pc C D {d D - a D ) C D d D (1 - V D ) 1 } 

where rj = d is the aspect ratio, and Co id the D-dimensional volume of a 
unit sphere. Hence, in mean field theory, p c is a monotonic function of the 
aspect ratio rj. This contradicts the result of simulations which show p c to 
have a pronounced minimum in 3 dimensions. However, some computer sim- 
ulations |9j show that as the dimensionality increases, the minimum becomes 
weaker, and at high dimensions, it appears that p c does increase ever more 
monotonically. Hence the mean field result corresponds again to the limit 
D -> oo. 

The importance of mean field theory is thus not so much in the numerical 
results is yields, but rather in the proof of the usefulness of the theoretical 
approach presented here. The mean field equations have been derived wholly 
within the Potts fluid picture, where heavy use is made of the available Hamil- 
tonian formulation. The s — > 1 limit then yields the percolation quantities, 
which do indeed have the expected mean field properties. It is hard to see 
how such results could have been obtained directly in the percolation picture. 
This indicates that other methods which can be applied to the Potts fluid 
should yield interesting results for the continuum percolation problem. This 
will be the subject of future inquiries. 



Appendix 

The functional F t [B] is defined as in Eq. (|3.8| ) 

F t (B) = IlnZ + ^^(g(z,i)exp[-^(i,i)]-l) 



22 



+ N_ 1) ^ ^ ^ [-f3v{i, j)\ 5 Xi ^ ) q 
+ ({3[h(j)-B(j)}^(\ 3 )} (A.l) 
We now have that 

In Z = N ln(r) - ln(AH) (A.2) 

where 

T = J2 [drexp [(3B(r)ijj(a)] (A.3) 
Similarly, we have immediately that 

(P IHj) - B(j)} ^(A,)) = fEfdj eP*™**) [h(j) - B(j)} (A.4) 

And 

(q(ij)exp [-Pv(i,j)]) 
= f^E j did 3 g(2,j>-^\^ B(i) ^ )+B( ^' )] (A.5) 

j)exp[-/3v(i,j)]5 A . )A .) o 
= ^ E J dido p(i, J >-^V^^) +fl W^)] ( 5 AilAj (A.6) 

Taking now the functional derivatives of these terms, we find 

Ml ' Z|) Lxp [/3£(£)V(A)] (A.7) 



5[iV/?.B(£)^(A)] r 

^[fe(j)-B(j)]V>(A 3 -)) 
5[/35(£)V(A)] 

= e^^wi {-1 + [h(x) - B(x)} ^(A)} 



e /3B(S)V(A) 



^ E / e^»<*'> [h(j) - B(j)\ ^(A,) (A.8) 



23 



5 (q(ij)exp [-(3v(i,j)] - 1)„ 

-Pv(x-y)3[B(y)i,(\ t ) 



5[(3B(x)^(X)\ 



A* 

2 



- e m^wA_Y, J didj q (iJ) e -^ e^WM+m^)] (A9) 
5 (p(i,j) exp [-(3v(i,j)] Sx^Jq 

(A.10) 

We now expand these expressions around B — > 0. We have that 

r = [drexp \j3B{r)il){ot)} « £ / dr [1 + /?S(f)^(«)] (A.ll) 

We note that 

2>(a) = ( a -l) + ( a -l)(-l) = (A.12) 
Therefore, to first order, 

r = sV + 0(B 2 ) (A.13) 

The identity Eq. (|A.12|) simplifies greatly all the functional derivatives. To 
first order, 

d ln Z " 4? ex P [^(^^(A)] + 0(5 2 )(A.14) 



5[N(3B{x)ijj{\)] sV 
5 (q(i,j)exp [-/3v(i,j)] - 1) 



5[/3B(x)i/j(\) 

And 

8[(3B(x)iP(\)] 



0(B 2 ) (A.15) 
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= e^^W-L {_1 + [h(x) - B(x)\ ip(X)} + 0(B 2 ) 

o V 

5 (p(i, j) exp [-/3v(i, j)} 5 Xi , Xj ) Q 

= ePB(£WX) -^T2 j dyp{x-y)e-^ s -®B{y) + 0(B 2 ) 

Summing up all these expressions and equating 5F t /5[(3Bip] 
obtain the condition 

B{x) = h(x) + ?-J dyp(x-$)e-^-®B(y) + 0(B 2 ) 
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